Charge transport and shot noise in ballistic graphene sheet 
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The current and the shot noise in a graphene sheet were analyzed in the ballistic regime for 
arbitrary voltage drops between leads and the sheet in the limit of infinite aspect ratio of the sheet 
width W to its length L, when quantization of transversal wave vectors is not essential. The cases 
of coherent and incoherent ballistic transport were compared. At high voltages the difference with 
coherent transport is not essential. But at low voltages conductance and Fano-factor dependences 
for incoherent transport become non-monotonous so that the conductance has a minimum and the 
Fano factor has a maximum at non-zero voltage bias. 
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Transport and shot noise are in the focus of recent experimental and theoretical investigations of graphene since they 
are tightly connected with unique spectrum of excitations in this material: the two-dimensional relativistic spectrum 
with the Dirac point. The presence of the Dirac point was experimentally confirmed in graphite (multilayered 
graphene )i and graphene^. From the very first experiments on graphene^ it was revealed that graphene has the 
minimum of conductivity, which is on the order of conductance quantum e 2 /h (see also Ref. Hf). Katsnelson^ and 
Tworzydlo et have shown theoretically that that the conductivity minimum must appear even in the ballistic 

S regime, i.e., in the absence of disorder. Tworzydlo et al£ have also analyzed shot noise and shown that close to the 
Dirac point the Fano factor of the graphene sheet coincides with that of the diffusive wire (F = 1/3). This stimulated 
'"^J i experimental investigations of shot noised. 

The analytical results of Katsnelsor£ and Tworzydlo et aif± were obtained for the case when the Fermi levels 
of leads are very close to the Dirac point of the sheet, and the main contribution to current and shot noise came 
i from evanescent states, i.e., from states classically forbidden in the sheet. In the case of the large voltage drops 

between the leads and the sheet the main contribution to the transport current and the noise power comes from 
t-H ■ the propagating modes, and the numerical calculations by Tworzydlo et al& for the finite aspect ratio W/L = 5 
demonstrated oscillations of the current and the noise power as functions of the large voltages. The question arises, 
whether these oscillations are connected with the finite aspect ratio, i.e. whether they remain in the limit W/L — » oo. 
Also, since the oscillations are related with coherent transport across the graphene sheet, they must be absent if 
transport is incoherent, i.e., a particle propagating in the sheet does not keep information on its quantum-mechanical 
phase. The present Communication focuses on these two issues: current and shot noise at large voltages for the 
infinite aspect ratio W/L, and comparison of coherent and incoherent ballistic transport. 

The analysis has been done for the case when leads are also made from graphene but strongly doped so that they 
have large Fermi wave vectors far from the Dirac point. Katsnelson^ and Tworzydlo et al£ considered exactly this 
case, and we can use their boundary conditions at two contacts (continuity of two-component electron wave function). 
Contacts are supposed to be ideal, and the only source of reflection is a mismatch of normal group velocities in the 
lead and in the sheet. Recent detailed theoretical investigations of contacts of graphene with various types of leads&i£ 
Vh ' have demonstrated insensitivity of transport through graphene to detailed properties of leads. 

We consider the graphene sheet of the length L and the width W with the large aspect ratio W L. The Fermi 
levels of the left and the right lead are shifted with respect to the Dirac point in the graphene sheet with the energies 
eVi and eVi respectively, so the voltage bias, which drives the current, is V = V\ — V% (see Fig. [J). The Hamiltonian 
of the graphene is well known: H = fivp(a x p x + a y p y ), where Up is the Fermi velocity, a x and a y are the Pauli 
matrices of the isospin, and p x and p y are inplane components of the momentum operators. The solutions of the 
Dirac equation are the two-component spinors, and the wave function for the electron crossing the left contact from 
the sheet (x > 0) to the left lead (x < 0) is 
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H ( _ K }+ ikn ) e -«--HM a t a, < 0. (1) 

Because of strong doping in the leads the values of the wave vectors in the leads K k, K x essentially exceeds the 
components k x and k y of the wave vectors in the sheet. On the other hand, because of the translational invariance 
along the axis y K y — k y . The reflection and transmission amplitudes are determined from the boundary conditions 6 
of continuity of the two spinor components. In the limit K k they are given by 

kx k iky ^ 2\fk x k . , 

k x H - k ~ j— iky kx k — t— iky 

and do not depend on K. This is another evidence that ballistic transport in graphcne is insensitive to details of 
leads. 

Similarly to this we can find the wave function for the electron crossing the right contact from the sheet (x < L) to 
the right lead (x > L). The reflection and transmission amplitudes in this case are r 2 = r^e 2lfcxi and t 2 = t\e % ^ kx ~ Kx ^ L . 
According to the theory of a double contact the total transmission amplitude for phase-coherent transport from the 
left to the second lead is 

ti*2 4k x k 
1 - rir 2 e H - i (k + k x ) 2 + fcg - [(k - k x ) 2 + k%]e 2ik * L ' V ; 

whereas the probabilities of transmission (T — \t\ 2 ) and reflection (R = \r\ 2 ) are 

\ht 2 \ 2 = H 

|1 - ri r 2 e 2ik ' L \ 2 kl + k 2 sin z {k x L) 
The analytic continuation of this expression yields the transmission probability for the evanescent mode k x = ip 

■2 



t=i-r= u ru a = „ , (4) 



j 1 E . (5) 

p 2 + k 2 sinh 2 {pL) 

These expression are known from Refs. Handi. The charge transfer via evanescent modes is in fact a direct single 
quantum- mechanical tunneling through the the whole sheet, which is classically inaccessible because of small energy 
of the electron. 

The total current consists from the current through the upper conductivity band (electrons) and the current 
through the lower valence band (holes), i.e. above and below the Dirac point respectively. We consider the case 
V = V\ > 0, Vq, = 0, when the whole current flows through the upper band. First let us restrict ourselves with 
propagating modes. Their contributions to the current and the differential conductance are given by 
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2eW r eV / hv p rV W6«f) 2 -^ 2evpW p eV / nv " f k 



dkx / dk„VxT = — / dk dk,, 

v x n 2 J J v kl + k 2 v sin 2 (k x L) 

ev F W f eV/hVF k x dk x V + ^V 2 ~ (v F Tik x /e) 2 cos{k x L) 



o cos(kxL) V - y/V 2 - (v F hkx/e) 2 cos(k x L) 

djp = q L r eV / hv F dk x 

" ! " ' ^/(eV/hv F ) 2 - k 2 x [(eV/Tiv F ) 2 sm 2 (k x L) + k 2 x cos 2 {k x L)] ' 
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Here Go = ^e 2 v F W /irhL is the minimum conductance, and v x = de/hdk x — v F k x /k is the x-component of the group 
velocity determined by the linear spectrum e = hv F k. Because of large aspect ratio W/L we replaced the sum over 
transversal components k y by an integral. In this limit the result does not depend on boundary conditions on lateral 
sheet edges discussed in Refs. [H and[H- 



At high voltages V 3> hv F /eL 
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At low voltages the transport via evanescent modes is more important. The expressions for the current and the 
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conductance through evanescent modes are (p 2 — k 2 — k 2 ) 
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2eW f eV / h -»F f°° 2ev F W f eV / h ^ f°° 
j e = —5- I dk x I dkyV x T = 5 — / dk dk y 



71 Jo J \/(eV/hv F ) 2 -kl 71 Jo Jk P 2 + k 2 . sinh (pL) 

ev F W f°° pdp y/(hv F p/e) 2 + V 2 cosh(pL) + V 



o cosh(pi) y/(hv F p/e) 2 + V 2 cosh(pi) - V 

— = GnL [°° p3dp (10) 

dV ° J ^(eV/hv F ) 2 + p 2 [{eV/hv F ) 2 sinh 2 { P L)+p 2 cosh 2 (pL)]' 



At high voltages V — > 00 this current saturates at the value 



2ev F W f°° pdp , , pL 1.82ev F W . , 

^ = —^1 ^ Z ylncoth T « ^ . (11) 
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Now let us write down the expressions for shot noise power. The contribution from evanescent modes: 

Ae 2 W f eV / hv " m f 00 f> m/i ms 4e 2 w F W^ /-^/^^ fj /- 00 ^ p 2 *; 2 sinh 2 (pi) 
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S e = 5— / cfe / dk v v x T(l — T) = £ — / djfc / dfc„ y ^— — . (12) 

' J ^eV/nv F )'-kl ^ 2 7o A y b 2 + A;2sinh 2 (pL)]2 V ; 



The propagating modes contribute the power 



ie 2 W f eV/nvF . u rVWJn^kJ ^ VpW r*v/*»* f k k 2 x k 2 y sin 2 (k x L) , 

Sp= 5— / dfe x / dk y v x T{\ — T) = / dk / dk y — ° , (13) 

7T 2 7o Jo Jo Jo y [k 2 x + k 2 y sm 2 (k x L)] 2 

At high voltages V — » 00 



Solid lines in Fig. [2] show the total conductance G = d(j p + j e )/dV and the Fano factor F = (S p + S e )/2e(j p + j e ) as 
functions of V = V\. The Fano factor changes from 1/3 at low voltages to 0.125 at high voltages. 

As well as in Refs. | and l6|, we assumed up to now that a particle propagated along the sheet coherently without 
any dephasing. Meanwhile, even in the ballistic regime the electron propagating across the sheet not necessarily 
"remembers" its phase. It is knownii that in the absence of the phase coherence the total transmission through the 
double contact is combined not from amplitudes of transmission but from probabilities of transmission. Namely, for 
the incoherent case one should replace Eq. ([5]) by 

T=l-R= TlT2 (15) 
1 - R1R2 k' v ' 

where T1.2 = | ^ 1 ,2 1 2 and R\ t 2 — | ^x,2 1 2 - Then the expressions Eqs. ([6]) and (fT"3"|) for the current and the shot-noise 
power should be replaced by 
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i* = —5r~ dk dk vir = - ^ It- , (16) 



k 7T 2 \Hv F J 



„ 4e 2 v F W [ eV ' hVF f k ,k x (k-k x ) 0.22e 2 v F W ( eV V 

Sp = Jo dk J dky k 2 = — [n^J ■ (17) 

The conductivity and the Fano factor for the incoherent propagation are shown in Fig. [2] with dashed lines. The 
comparison of two cases shows that at high voltages the role of possible phase coherence is not essential, and the 
oscillations resulting from the phase dependence are insignificant as was observed in the experiments^. Apparently 
the reason of it is that even in the phase-correlated case there is decoherence from summation over the states with 
different wave vectors. The most remarkable difference between the two cases is non-monotonicity of dependences in 
the incoherent case at low voltages: the conductiance has a minimum, where the conductance is slightly less than 
the conductance at V = 0, (so the latter is not the minimum conductance anymore!), whereas the Fano factor has a 
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maximum, where its value 0.37 exceeds 1/3. It is worthwhile to note that DiCarlo el alZ observed the Fano factor 
up to 0.38, which could be an evidence of incoherent propagation in the sheet. 

A similar contribution can be done for nonzero voltage V = V 2 (with V\ = 0) independently on whether the charge 
is transported only by electrons (V 2 > 0) or also by holes (V 2 < 0). Symmetry between electron and holes is valid as 
far as the Fermi wave vectors in the leads essentially exceed the wave vectors k inside the sheet. In general the total 
current through the sheet and the Fano factor are determined by the relations 

J = j P (Vi) +UVi) -j P (V 2 ) - Je (V 2 ) (18) 

and 

1 SpM) + SeW - S P (V 2 ) - S e (V 2 ) 

2e jp(Vx) + j e (Vi) - j p (V 2 ) - j e (V 2 ) ' 1 ' 

One might expect that the ballistic transport itself should be lossless and noiseless. But it is well known that 
transport through any ballistic conductor is inevitably accompanied by dissipation and noise at the contacts^: elec- 
trons, which are injected from one lead to another, cross the ballistic conductor without scattering, but disturb the 
equilibrium in the lead, where they are injected. This results in dissipation inside leads. Resistance and noise for the 
transport through graphene also originate from stochastic processes of relaxation in the leads near the contacts. But 
the unique feature of graphene is that the contact phenomena are governed by the peculiar spectrum in the Dirac 
point rather than by details of the spectrum in leads. 

In summary, we analyzed the current and the Fano factor as functions of voltage in a graphene sheet for arbitrary 
voltage drops between leads and the sheet in the limit of infinite aspect ratio of the sheet width W to its length L. 
The oscillations on these dependences, which were revealed earlier 6 for the finite aspect ratio, are insignificant in the 
limit W/L — > 00. The case of incoherent ballistic transport was also analyzed. At high voltages the difference with 
coherent transport is not essential. But at low voltages conductance and and Fano-factor dependences for incoherent 
transport become non-monotonous so that the conductance has a minimum and the Fano factor has a maximum at 
non-zero voltage bias. 
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by the Large Scale Installation Program ULTI-3 of the European Union. 
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FIG. 1: (color online) Energy bands in the left and the right leads and the graphene sheet (center), a) Voltage bias Vi > is 
applied between the left lead and the Dirac point in the sheet, Vi = 0. b) Voltage bias V2 < is applied between the Dirac 
point in the sheet and the left lead, V\ = 0. c) Both the voltage drops, Vi and V2, are present. The darker areas in the sheet 
energy pictures show the fully occupied states in the sheet. The lighter-shaded areas show the partially occupied states (only 
left-moving particles and right-moving holes). 
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FIG. 2: (color online) The relative conductance G/Go (a) and the Fano factor F (b) as functions of V — Vi (V2 = 0). Solid 
lines refer to the phase-coherent propagation through the sheet, and dashed lines refer to the phase-incoherent case. 



